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Abstract: In a recent work, S. Dowker has shed doubt on a recipe used in computing the par- 
tition function for a matrix valued operator. This recipe, advocated by Benson, Bernstein 
and Dodelson, leads naturally to the so called multiplicative anomaly for the zeta-function 
regularized functional determinants. In this letter we present arguments in favour of the 
mentioned prescription, showing that it is the valid one in calculations involving the rela- 
tivistic charged bosonic ideal gas in the framework of functional analysis. 



PACS numbers: 05.30.-d,05.30.Jp,11.10.Wx,11.15.Ex 

In a recent work, Stuart Dowker Q has shed doubt on a commonly used manipulation in 
computing functional determinants related to matrix valued elliptic operators. This recipe is 
widely used 0-0] and can be summarized as follows. Within the one-loop approximation or in 
the external field approximation, one often has to evaluate Euclidean functional integrals of the 
kind 

Z = J DfaDfoe' I , (1) 
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where A is a matrix valued differential operator, which we assume to have constant coefficients 
(this is not a restriction as long as one has to compute the one- loop effective potential). The 
result of the Gaussian functional integration is 

Z = (det Ay 1/2 . (2) 

The problem is how to compute the functional determinant in Eq. (|2|). Note that A has matrix 
elements with discrete (field) and continuous (spacetime) indices. If A is also diagonal in the 
discrete indices, one has 

In his work Dowker presents two different ways of computing this determinant, claiming 
that one of them does not lead to correct results. The first recipe || is equivalent to taking the 
algebraic determinant first and the functional one afterwards, i.e. 

(detA)i = det(LiL 2 ) . (4) 

According to Dowker jC], it could seem more natural to use the alternative recipe 

(det A) 2 = detLidetL 2 . (5) 

This is seen as the implementation of the right way to take the functional determinant of the 
matrix valued operator, i.e. an ordinary determinant on both the spatial and field indices, at 
the same moment. The example analized is that of two real free scalar fields of different masses, 
for which the partition function has the form in (|l|) and the operator is diagonal in the field 
indices i,j as in (||). 

This second recipe seems quite natural indeed, since it is well known that for finite matrices, 
det^4-B = det A det B. Unfortunately, in the continuum, one needs a regularization and for 
regularized functional determinants the two recipes may give different answers, because of the 
existence of the so called multiplicative anomaly, discovered by Wodzicki (see, for example, ||), 

det AB = det A det Be a{A ' B) , (6) 

and, also in very simple cases of physical interest, it is possible to show that a(A, B) is not 
vanishing ||. Thus the question posed by Dowker is substantial. 

In our opinion, both recipes are formally acceptable. In fact if one considers the eigenvalue 
problem for the matrix valued operator A, one arrives at the formal determinant: det A = 
n nin2 A ni A„ 2 , which can as well be rewritten as det A = ]J A ni ]J n2 A„ 2 . Of course, as rigorous 
equalities these expressions are restricted to the finite dimensional case (and to a limited class 
of absolutely convergent situations). 

In the following we would like to present arguments in favour of the first recipe, analysing 
the general validity of both the recipes and also considering as a crucial example the case of a 
free relativistic charged bosonic field at finite temperature ||, [l0| . The self-interacting charged 
scalar field was studied in || . 

We have argued that both recipes actually coincide in the finite dimensional case. One 
should keep in mind that any extension to the continuous, functional case always starts from 
a discretization. In particular, the finite dimensional example posed by Dowker at the end of 
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]|] pretending to prove a discrepancy already at this level does not apply. In fact, let us start 
considering Dowker statements in general. He supports the inequivalence of the two recipes 
analysing a generic four by four matrix A a pij diagonal on where the indices a,j3 = 1,2 
represent the discretized version of the continuous space-time indices and i, j = 1,2 the fields 
indices. 

Then, the two recipes correspond to det a ^ detij A a ^ and deto-^^j A a ^ respectively. For 
such a generic matrix it is straightforward to see that the two results are different and Dowker's 
statement that the first recipe could give inexact results seems correct. The point, though, is 
that the matrices we encounter in field theory have additional structural requirements. In the 
continuous limit, the matrix valued differential operator is normally diagonal in the continuous 
indices, namely 



A i,j{x,y) = Aij5(x,y) . 



(7) 



and will therefore have a block-diagonal structure. The determinant of such a matrix is the 
product of the determinants of the blocks, therefore 



det A 



det det A, 

a,/3 i,j 



(8) 



which is exactly the widely used recipe. 

A more formal analysis would require the study of the proper discrete matrix, for the func- 
tional integral is solely defined as the continuum limit of a discretized lattice version. The 
derivative there is in fact defined as difference of the values of the field in two neighbouring 
points of the lattice and the derivative squared has, therefore, terms which are off-diagonal in 
the spacetime indices, like 4>(x + l)(p(x). In one dimension the corresponding matrix A a pij would 
have the structure 
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(9) 



where □ and O represent blocks in the field indices. It is easy to see how, for a matrix with such 
a structure, the above equality (||) holds again, since the only contributions to the determinant 
come from the diagonal blocks. 

It thus seems reasonable that the first, widely used recipe for calculating the functional 
determinant is equivalent to the other more rigorous one in the finite limit of field theory. This 
shows also that the reasons for the presence of the anomaly have to be found in the infinite 
nature of the functional determinant. 

Let us turn our attention to the free charged bosonic field at finite temperature. In order 
to evaluate functional determinants we will make use of zeta-function regularization [11, O], 



Recall that the one-loop Euclidean partition function, regularized by means of zeta-function 
techniques, reads IT! 



InZ 



- — In det 
2 



L D £< 



ic'(0|L D )-ic(0|L D )ln^ 



where (,(s\Ld) is the zeta function corresponding to Lp, a second order elliptic differential 
operator, £'(0\Lei) its derivative with respect to s, and t 2 is a renormalization scale parameter. 
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The grand canonical partition function for an ideal charged relativistic boson gas may be 
written as (// is the chemical potential) 



D(f>ie~^ fo ' dT I d3xd3 y<t ) i( x ) A iA x >y)<t > j(y) ) ng) 

(t)=0(t+/?) 



where the two real fields (pi, i = 1, 2, are chosen to describe the degrees of freedom of the boson 
gas, and the operator A has matrix elements given by 

Aij(x,y) = [hij + 2/j,e ij y r L?J S(x,y), (11) 

with 

L fi = (L r + L 3 - /i 2 ) <%, L 3 = -A 3 + m 2 , (12) 

in which A3 is the Laplace operator on IR 3 (continuous spectrum k 2 ) and L T = — S 2 (discrete 
spectrum u 
written as 



spectrum w 2 = 47 ^" ) the Laplace operator on S 1 . In this case, the partition function may be 



The first recipe consists in taking first the algebraic determinant and then the functional deter- 
minant. The result is 

Zfo = (det {i 4 [{L T + L 3 - /i 2 ) 2 + 4/x 2 L r ] })~ V2 = (det(^ 4 L+L_)) ^ , (14) 

where 

L± = L T + L 3 + /i 2 ± 2/i (L 3 )3 . (15) 

In an attempt of using the second recipe, one may observe that every 2x2 matrix here can 
always be diagonalized. Then, modulo a trivial functional Jacobian corresponding to the diag- 
onalization, one has 

and the answer coming from the second recipe is 

Z^ = (det(£ 2 L + )det(£ 2 L_))~ 1/2 . (17) 

On the other hand — as is usually done in quantum field theory — one can also describe the 
gas by two complex fields <j> and <p* , defined by 

= -L(0 1 + »^) j>* = -L^-ifc) . (is) 

The corresponding grand partition function reads Q 

Z» -(**{*>(«+ 1)}Y' 2 (19) 



-1/2 
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where now 



Note that we have L + L 
approaches, namely 



K± = L 3 + L T - i? ± 2»/z (L T ) 2 . (20) 
. = K+K-, thus the first recipe gives the same answer for the two 



Zfo = (det [l\L+L_ 
The second recipe gives, on its turn, 

Z p>ft = (det det (l 2 L_ 



-1/2 



-1/2 



(det [i 4 (K + K_) 



-1/2 



(det (£ 2 K + )det {fK-)) 



-1/2 



(21) 



(22) 



When the multiplicative anomaly is non vanishing, one of two recipes is in contradiction. In Ref. 
[13] it has been shown that in odd dimensions, the multiplicative anomaly is vanishing and the 
two recipes give the same answer. In even dimensions, in particular in H 4 , the multiplicative 
anomaly is non vanishing and only the first recipe gives the same partition function for the two 
approaches, since we have 

In (^ 2 detL + ) + In (^ 2 detL_) + a(L + , L_) = In (l 2 det 2f+) + In (fdetK^ + a(K + , KJ) (23) 
but, on the contrary, 



In (l 2 det L+) + In (t 2 det L_) / In (f 2 det #+) + In (l 2 det K. 



(24) 



It has also been shown that the statistical sum contribution to the grand partition function is 

ression [10] 

-/3(v%+a0) + l n (1 _ p-flv^-^ 



the same and yields the well known expression [10] 



s(/M = £[i*(i 



i 



(25) 



where = k 2 + m 2 . The discrepancy manifests itself in the " vacuum sector" , namely in the 
contribution to the grand partition function linear in (5 and the presence of the multiplica- 
tive anomaly, first recipe, renders the grand partition functions equal and independent by the 
parametrization of the degrees of freedom of the relativistic ideal gas. 

Having performed the above calculation — that seems to leave little chance for discrepancy — 
one might still ask: how could the apparently more natural second recipe Eq. (|5|) fail ? Well, 
the answer is, to begin with, that we are dealing with a very elusive mathematical and physical 
point. In fact, for a direct sum of operators, acting on a corresponding direct sum of independent 
functional spaces (this is the case in quantum physics when, for instance, a superselection rule 
is imposed upon the system), we indeed have the factorization property: 



det A 



( Al 


V i 





A 2 



det Ai det A 2 



(26) 



In this situation the product itself of the operators, AiA 2 , makes no sense in general, much less 
its determinant, and the prescription det Ai det A 2 • • • is to be used. But things are absolutely 
different when one is working within a functional space where field mixing and rotations are 



5 



allowed, and one just obtains the diagonal expression for the matrix of operators acting in this 
space as a particular form, after a convenient diagonalization process (this is precisely what 
happens in our example above). The moral we extract from the outcome of our analysis is 
that the invariant which is preserved under this process of change of basis is precisely the 
determinant of the operator matrix, but not the fact that it is equal to the product of the 
functional determinants of the operators. It is precisely that invariance what lays in the heart 
of our example: the noncommutative anomaly is the missing term necessary in order to preserve 
it. 

To summarize, we have here carried out what is, in our opinion, a serious consistency check 
in favour of the first recipe for the calculation of determinants of matrix valued operators and, as 
a consequence, provided arguments in favour of the relevance of zeta-function regularization and 
related multiplicative anomaly in quantum field theory. As far as this last issue is concerned, in 
Ref. []14] we also respond to a criticism appeared recently [T^] , concerning the use of zeta-function 
regularization. 
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